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INVESTIGATION  OF  THE  PROPAGATION  OF  BLAST  WAVES  OVER  RELATIVELY 
1  LARGS  DISTANCES  AND  THE  DAMAGING  POSSIBILITIES  OF  SUCH  PROPAGATION 

V  ABSTRACT 

The  t^st  firing  of  large  guns  and  explosives  has  resulted,  on  numerous  occasions,  in  complaints 
from  inhabitants  residing  outside  the  testing  area.  Tfie  size  of  the  testing  ar.ea  should  preclude  any  such 
damaging  effects  but  the  presence  of  these  effects  indicated  that  on  certain,  days  a  blast-wave  is  propa¬ 
gated  to  unusually  great  distances  witha  relatively  high  pressure  difference  across  the  wave  surface.  Con¬ 
sideration  of  seismic  waves  and  air  waves  lgd  to  the  conclusion  that  this  unusual  blast  effect  iq,  propagated 
through  the  atmosphere  and  hence  is,  to  a  great  extent,  dependent  upon  existent  meteorological  conditions. 

In  this  paper  an  attempt  is  made  to  determine  those  meteorologi  al  conditions  which  are  conducive 
to  unusual  blast  wave  propagation  and  to  establish  a  few  rules  of  thumb  which  will  enable  a  reasoning - 
son,  with  the  aid  of  temperature  and  wind  soundings: of  the  lower  atmosphere,  to  determine  the  presence  of 
dangerous  propagating  conditions.  These  rules  cf  thumb  are  not  based  upon  a  rigid  theoretical  analysis 
simpiy  because  our  lack  of  knowledge  of  microscopic  meteorological  phenomena  and  the  physical  charac¬ 
teristics  of  spherical  blast  waves  would  make  such  a  rigid  analysis  meaningless. 

The  work  of  various  individuals  in  the  field  of  sound  propagation  through  the\tmosphere  is  stuaied 
and  the  results  of  their  work  are  employed  in  the  analysis  of  the  subject.problem.  while  initially  a  blast 
wave  may  differ  greatly  from  an  ordinary  souhd  wave,  a  short  distance  from  the  explosiVesource.these  two 
types  of  wavesttravel  with  identical  velocities  and  are  affected  in  equal  degrees  by  exismng  meteorological 
conditions.  \ 

Because  of  mathematical  and  physical  difficulties  the  necessary  approach  to  the  problem  was  to 
study  certain  Ideal  atmospheric  structures  and  to  determine  the  altitude  function  of  temperature  and  wind 
(and  relative  humidity)  which  is  necessary  if  destructive  blast  propagation  is  to  take  place.  The  results  of 
these  studies  indicated  that  a  necessary,  though  not.-sufficient  condition  is  that  a  temperature  inversion 
exists  in  the  atmosphere  through  which  the  blast  wave  is  propagated;  i.e.,  an  increase  in  temperature  -with 
altitude; 

The  necessary  increase  in  temperature  with  an  increase  in  altitude  to  produce,  destructive  blast 
propagation  immediately  suggests  certain  meteorologies  phenomena  associated  with  this  temperature  con¬ 
dition.  'Employing  these  meteorologic  phenomena  enables  one  to  set  forth  certain  rules  of  thumb  which  will 
indicate  dangerous  blast  conditions. 

1.  A  high  pressure  area  (atmospheric)  which  is  stationary  over  the  area  in  question  becomes  danger¬ 
ous  after  the  first  day  of  its  presence. 

2.  Large  diurnal  temperature  variations  at  the  earth’s  surface  indicate  unfavorable  firing  conditions.. 

3.  Poor  visibilities  and  light  winds  in  the  early  morning  followed  by  cloudless  conditions  through¬ 
out  the  day  indicate  dangerous  blast  propagation. 

4.  Light  surface  winds  and  low  relative  humidity  at  the  ground  plus  the  presence  of  stratified  cloud 
deck  below  10,000  feet  indicate  dangerous  firing  conditions. 

Poor  visibilities  and  light  winds  in  the  early  morning  followed  by  the  formation  of  cumulus  (swell¬ 
ing)  clouds  in  the  late  morning  and  early  afternoon  indicate  safe  conditions  for  an  afternoon  firing. 

Several  situations  in  which  destructive  blast  was  reported  were  picked  at  random  and  the  meteoro¬ 
logical  conditions  occurring  at  that  time  were  examined.  It  was  found  in  each  of  these  cases  that  a  positive 
vertical  temperature  gradient  occurred  in  the  lower  atmosphere  (below  10,000  feet)  and  that  the  weather 
conditions  fit  one  or  more  of  the  4  situations  described  above. 
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INTRODUCTION 

Ac  Integral  portion  of  the  research  carried  on  at  the  Aberdeen  Proving  Ground  Involves  the  study 
of  statically  detonated  charges.  Since  some  of  these  charges  must  of  necessity  be  quite  large,  the  blast 
waves  produced  by  the  explosions  have  large  pressure  discontinuities  across  their  boundrles  and  are  cap¬ 
able  of  inflicting  serious  damage  upon  objects  a  short  distance  fromthe  pressure  wave  source.  However, 
beyond  a  relatively  short  distance,  the  Intensity  of  the  pressure  or  blast  wave  decreases  approximately  as 
the  square  of  the  distance  from  the  source,  and  the  area  cf  the  testing  range  Is  such  that,  .theoretically,  no 
blast  damage  should  occur  to  objects  without  the  confines  of  this  area.  Unfortunately,  blast  damage  has 
qccured  in  the.  past  —  at  considerable  expense  to  the  Government  —  and  It  seems  highly  desirable  to  In¬ 
vestigate  the  conditions  under  which  this  destructive  blast  Is  propagated.  During  the  war  the  press  and 
importance  of  blast  Investigations  —  of  charges  dropped  from  aircraft  as  well  as  static  detonations  —  for¬ 
bade  consideration  of-  minor  damage  to  private  property  and  the  policy  of  firing  first  and  paying  afterwards 
was.  generally  adopted.  With  the  cessation  of -hostilities  came  budget  reductions  and  It  became  at  once  de¬ 
sirable  to  find  some. method  of  predicting,  beforehand,  the  presence  of  unfavorable  meteorological  con¬ 
ditions,,  i.e.,  conditions  when  blast  damage  Is  likely  to  occur. 

The  phenomenon  of  unusual  sound  or  blast  propagation  has,  of  course,  been,  known  for  many,  years. 
It  has  been  the  subject  of  many  Investigations  during  the  past  century  by  such  men  as  Stokes  (1),  Reynolds 
(2),  Milne  (3),  and  Whipple  (4)  In  England  and  Gutenberg.  (5)  and  Duckert  (6)  In  Germany.  Successive  zones 
of  audibility  and  silence  occuring  when  severe  explosions  have  taken.place  have  been  plotted  and  analyzed 
and  their  causes  determined;  During  the  first  World  War' the  sound. ranging  of  artillery  fire  was  used  for 
determining  enemy  positions  and  it  was  quickly  noted  that  on  occasions  the  intensity  of  recorded  gunfire 
was  quite  great. 

During  the  months  of  November  and  December  of  the  year  1940,  the  US  Bureau  of  Mines  (7)  con¬ 
ducted  a  'series  of  experiments  at  the  Aberdeen  Proving  Ground  In. which  seismic  and  air  pressure  meas¬ 
urements  were  made  in  coincidence  with  the  firing  of  relatively  low  calibre  guns.  For  the  experiments 
several  houses  of  moderate  size  and  construction  located  on  the  northwesterly  boundry  of  the  Proving 
Ground  were  chosen.  Seismometers  were  mounted  on  the  floors  of  these  dwellings  and  measurements 
taken  while  firing  was  In  progress.  Results  of  the  experiments  Indicated  that  the  firing  of  8  Inch  guns  and 
8  inch  Howitzers  produced  no  superficial  damage  to  these  houses  located  only  1  to  2  miles  distant  from 
the  point  of  fire.  No  further  experiments  along  these  lines  were  conducted  by  the  Bureau  of  Mines. 

Iri  June,  1944,  R.  G.  Sachs  (8)  of  the  Ballistic  Research  Laboratories  published  a  report  entitled, 
“The  Limiting  Distances  for  Superficial  Damage  to  Buildings  due  to  the  Firing  of  Guns  and  High  Ejq>losives." 
In  this  report  the  author  reviewed  the  old  Bureau  of  Mines’  experiments  in  the  light  of  the  complaints  then 
being  received  from  Inhabitants  residing  near  the  Proving  Ground.  A  survey  of  these  complaints  and 
claims  against  the  Government  for  damage  caused  by  firing  programs  enabled  Sachs  to  evolve  an  empirical 
formula  which  defined  the  limiting  distance  beyond  which  destructive  blast  from  detonated  charges  evident¬ 
ly  did  not  occur.  This  formula  is  given  as  o 

~>R  =  3750  VW 

9  ft  =  Llmltlng^lestructlve  range  (feet) 

(  W  =  Weight  of  explosive  (pounds) 

<,  =1/7  weight  of  propellant  charge  if  gun  is  used. 
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Re-examlnatLon. of  the  data  gathered  In. the  Bureau  of  Mines’  experiments  revealed  that  In -light  of 
the  above  formula  no  destructive  effects  from  the  gunfire  could  be  expected.  The  experiments  did  show,  __ 
however,  *hat  the  transmission  of  shock  through -the  ground  was  of  ho  more  importance  than  the  transmis¬ 
sion  of  blast  through  the  atmosphere. 

The  ever  present  danger  of  destructive  blast  propagation  from  detonated  charges  resulted,  quite 
naturally,,  in  the  formation  of  a  central  uthoritative  source  which  had  the  unpleasant  task  of  determining 
whether  or  not  a  firing  should  take  place  i.e.,  whether  or  not  destructive  blast  was  likely  to  occur.  This 
authoritative  source,  the  Safety  Officer,  could  only  base  his  judgments  on  analogous  situations  when  un¬ 
usual  blast  propagation  had  been  known  to  occur,  and,  because  the  responsibility  was  great  and  the  cri¬ 
teria  meager,  much  time  and  money  have  been  wasted  by  overcautious  delays  in  firing  programs.  Actual¬ 
ly  it  is  not  at  all  unreasonable  to  suppose  that  more  government  funds  have  been  expended  in  firing  delays 
than  would  have  been  used  had  the  firings  proceeded  as  scheduled  and  payments  made  for  any  reasonable 
claims  for  resulting  damage.  It  most  certainly,  is  a  fact  that  complaints  have  been  received  (although  no 
monetary  claims  were  involved)  for  a  firing  which  was  believed  "safe11. 

It  must  be  polntea  out  that  no  reflection  upon  the  ability  of  any  individual  is  implied  in  these  state¬ 
ments.  The  basis  of  prediction,  built  upon  the  weakest  of  foundations,  is,  of  course,  at  fault-. 

The  investigation  of  unusual  blast  propagation  through  the  atmosphere  is  detailed  in  the  following 
pages  and  some  attempts  are  made  to  establish  certain  fundamental  rules  for  producting  a  more  accurate 
prediction  of  dangerous  propagating  conditions, 

THEORY 

The  propagation  of  sound  through  the  atmosphere  is  a  well  understood  physical  phenomenon- mid  is 
quite  capable  of  being  mathematically  analyzed.  The  wave  and  wave  front  characteristics  of  a  small  pres¬ 
sure  disturbance  immediately  suggest  a  Huygen’s  construction  for  determining  future  frontal  positions  of 
such  a  disturbance.  Vertical  or  horizontal  wind  and  temperature  gradients  which  exist  in  the  atmosphere 
may  be  simply  treated  by  considering  the  propagation  velocity  to  vary  along  the  wave  surface  of  a  small 
pressure  disturbance  and  succeeding  positions  of  the  wave  front-plotted  accordingly..  Actually,  in  many 
respects,  phenomena  in  sound  propagation  through,  the  atmosphere  are  quite  analogous  to  the  better  known 
and  better  understood  optical  counterparts.  Refraction,  reflection,  and  diffraction  are  common  occur¬ 
rences  in  either  type  of  propagation  and  their  causes  may  be  readily  determined  by  an  examination  of  pre¬ 
vailing  meteorological  conditions.  The  analogy  ceases,  however,  when  consideration  is  taken  of  the  fact 
that  density  differences  rather  than  temperature  variations  affect  the  passage  of  light  through  the  atmos¬ 
phere.  Thus  In  a  normal  atmospheric  temperature  structure  (temperature  decreasing  with  altitude)  light 
will  be  refracted  downward  toward  the  region  of  greater  density  whereas  sound  will  be  refracted  upward 
toward  the  region  of  lower  propagation  velocity.  Also  wind  gradients  (absolute  values)  have  a  negligible 
effect  upon  light  propagation  but  are  most  important  in  the  treatment  of  sound  wave  motions. 

Lord  Rayleigh  (9)  has  furnished  a  most  useful  tool  for  the  study  of  sound  propagation;  namely,  the 
concept  of  "rays  of  sound."  Although  the  concept  may  be  physically  artificial,  the  location  of  sound  fronts 
by  suitable  velocity  vectors  (sound  rays)  and  Huygen’s  constructions  is  most  accurate.  In  fact,  the  pro¬ 
pagation  of  sound  through  a  medium  in  which  the  velocity  vectors  change  continuously  -  but  not  according 
to  theoretical  or  empirical  laws  -  may  only  be  treated  accurately  by  just  such  a  concept.  Figures  la  and 


FIGURE  la.  Shows  the  situation  giving  rise  to  refraction  of  sound  by  wind  effects  alone.  Uj,  u2,  and 
u.  are  the  varying  wind  velocities.  C  is  the  velocity  of  sound  arising  from  a  uniform  temperature  in  the 
medium.  Sj,  S2,  S3  and  S^  are  sound  wave  frontal  positions  after  time  intervals  of  6  tj,  6  t2,  6 13  and 
5  t^.  Figure  la  shows  quite  obviously  that  the  sound  ray  is  the  resultant  of  the  vectors  of  u  i  t  and 
c.a  t  and  also  that  it  is  not  normal  to  the  wave  front  in  the  general  case.  In  the  atmosphere  where  i  t 
must  of  necessity  be  taken  small,  the  sound  ray  would  appear  as  a  continuous  curve. 


FIGURE  lb.  Shows  the  refraction  of  sound  by  a  temperature  gradient  alone.  In  this  case  there  is  no  wind. 
Since,  for  the  purpose  of  illustration,  6 1  had  to  be  made  finitely  large,  the  sound  rays  are  drawn  as  straight 
lines  and  the  sound  velocities  Cj,  c2,  Cj  and  must  be  taken  as  average  velocities  over  the  altitude  inter¬ 
val  indicated.  In  the  atmosphere,  of  course,  the  sound  rays  would  appear  as  continuous  curves  since  i  t 
could  be  made  sufficiently  small. 
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lb  Illustrate  the  relationship  between  sound  raj'  and  sound  front  under  different  pr.pagailng  ;  mditi  .r_ . 

Thus  far  the  discussion  has  been  concerned  only  with  the  propagation  of  sound  through  the  atm  ,c  - 
pnere  whereas  the  problem  itself  involves  the  propagation  of  blast  or  shock,  waves.  Btuaose  a  shock  wave 
is  characterized  by  a  nearly  instantaneous  pressure  rise  along  Its  frontal  surface  and  because  the  exresc 
pressure  behind  the  shock-surface  may  even  be  whole  multiples  of  the  atmospheric  pressure  (close  i.  the 
origin  of  the  shock  wave),  the  atmosphere  is  heated  irreversably  (entropy  increase)  Ly  the  passage  A  such 
a  wave.  For  these  reasons  —  and  also  because  of  the  existence  of  Individual  atmospheric  particle  m:t:_r~ 
due  to  the  detonation  —  the  blast  disturbance  does  not  proceed  outward  at  the  velocity  of  sound  tat  rather 
at  a  velocity  dependent  upon  a),  the  size  of  the  detonated  charge,  b)  the  distance  of  the  blast  disturbance  from 
the  activating. source  and  c)  the  existing  meteorological  conditions.  However,  since  the  tutal  energy  A  a 
blast  wave  must  remain  constant  as  the  wave  proceeds  from  the  source,  the  pressure  behind  the  shock 
surface  decreases -.and  eventually  reaches  a  value  comparable  to  that  found  in  ordinary  ssund  AavcS. 
Thereafter  the  blast  wave  is  propagated  at  the  local  velocity  of  sound. 

The  rapidity  with  which  the  pressure  decreases  (as  a  function  of- distance  from  the  source)  behind 
a.shock  surface  is  dependent,  quite  obviously,  upon  the  size  of  the  originating  disturbance  ar.d  seems  i; 
be  aest  determined  oy  empirical  methpds.  Taylor  (1C)  has  discussed  the  behavior  of  a  blast  wave  emanat¬ 
ing  from  the  detonation  of  a  500  lb.  charge  and  gives  tables  showing  the  variation  of  pressure  be  hind  the 
shock  wave  as  a  function  of  distance  from  the  source.  Extrapolation  of  Taylor’s  results  (Fig.  2)  shows 
that,  at  500  feet,  the  velocity  of  the  blast  wave  (from  a  500  lb.  bomb)  is  less  than  1.01  times  the  velocity 
of  sound  at  the  same  atmospheric  temperature.  Over  distances  of  two  miles  or  greater  it  would  be  ex¬ 
tremely  difficult  to  distinguish  time* differences  of  arrival  between  sound  and  blast  wave-  originating  fr:m 
a  common  point  at  a  common  time.  Even  for  charges  of  considerably  greater  magnitude,  propagation  vel:- 
city  of  the  blast  wave  may  be  assumed  equal  to  that  of  sound  where  distances  involved  are  larger  than  tv;-; 
or  three  miles.  This  assumption,  while  not  entirely  valid,  enables  one  to  apply  the  propagation  laws  ox 
sound  to  blast  wave  phenomenon  —  about  which  far  less  is  known  and  understood  (referring  specifically  t: 
the  spherical  case). 

One  further  difference  between  sound  and  bl  ;st  phenomena  must  be  indicated  at  this  point.  The  wave 
length  of  a  blast  wave  changes  as  a  function  of  distance  from  the  source  (at  least  during  the  relatively  early 
expansion)  and  the  attenuation  suffered  by  a  pressure  wave  in  the  atmosphere  seems  to  be  dependent,  to  a 
large  extent,  upon  tne  frequency  of  the  disturbance.  Unfortunately  very  little  is  known  about  this  change  in 
wave  length  mentioned  above  and  further  assumptions  must  be  made.  Since  sound  waves  do  not  vary  in  fre¬ 
quency  along  a  path  it  seems  logical  to  assume  that  when  the  pressure  amplitude  of  a  blast  wave  approaches 
that  found  In  a  sound  disturbance  the  wave  length  tends  to  become  constant  and  undergoes  no  change  beyond 
a  certain  limiting  distance.  If  this  assumption  is  made  it  then  follows  that  blast  waves  originating  from  dif¬ 
ferent  size  charges  will  have  different  final  frequencies  (othe:  things  being  equal)  and  will  be  affected  in 
different  ways  during  passage  through  the  atmosphere.  Data  on  sound  ranging  of  guns  during  the  first  World 
War  (11) indicated  that  the  wave  length  of  the  pressure  <vave  from  gun-blast  was  proportional  to  the  calibie 
of  the  gun  and,  for  a  15  Inch  gun  the  frequency,  at  a  great  distance  from  the  source,  was  found  to  be  4  or  5 
cycles  per  second.  (Wave  length  defined  as  time  length  of  positive  pressure  phase). 
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Consideration  of  the  facts  and  suppositions  indicated  above  leads  to  theiollowlng  conclusions: 

1)  Blast  waves  may  be  treated  as  ordinary  sound  waves,. with  regard  to  propagation  and  physical 
Characteristics,  beyond  a  relatively  short  distance  from  the  source  of-disturbance. 

2)  Propagation  laws  developed  for  sound  phenomena  do  not  consider  attenuation  as  a  function  of 
frequency  or  as  a  function  of  water  vapor  and  dust  particles  in  the  air.  Therefore,  a  disturbance  of  greater 
wave  length,  such  as  a  blast  wave,  may  be  expected  to  conform  more  closely  to  the  idealized  formulae  than 
one  abthe  lower  or  intermediate  audible  levels. 

fit-1916,  Fujiwhara  (12)  published  the  second  of  two  reports  in  which  he  described  an  investigation  of 
the  sound  ray  paths  emanating  from  the  Asama  volcano  eruptions  of  1909  to  1913.  fit  treating  this  problem 
Fujiwhara  was  forced  to  make  two  basic  and  simplifying  assumptions: 

a)  only  wind  and  air  temperature  affected  propagation  velocities  and 

b)  horizontal  gradierts  in  wind  and  temperature  were  zero  (i.e.  wind  and  temperature  varied  only 
with  height)..  Thus  if  x  and  y  are  rectangular  coordinates  of  the  earth's  surface  considered  as  a  plane,  the 
projected  curve  of  a  sound  ray  on  this  plane  is  given  by 


tx'I 


'  c  cosip+  u^j  (c^sin  iQ)  -  (u^-uq)  cos  <p~  (v^-vQ)sLn  rdz 


c[  fcQ/sln  lQ)  -  (Uj-uq)  costp-  (Vj-v^sin 


V-  I 


c  sLntp+  Vj  j  (CQ/sin  iQ)  -  (Uj-uq)  cos  tfr-  (v^-vjshujf  jdz 


o  ct  J(cQ/sin  iQ)  -  (ux-u0)  cos  V- 


V-  (v^-v^jsin 


where 


z  =  altitude  -  vertical  axis; 

c  =  velocity  of-sound  due  to  temperature  alone.  (Equation  A.7  Appendix  A) 
u,  v  =  horizontal  wind  components- in  x  and  y  directions  respectively, 
i  =  angle  between  z-axls  and  normal  to  the  wave  front. 

=  azimuth  of  wave  front  normal  for  sound  ray  in  quest'on  with  respect  to  x-axis. 
Subscripts  0  and  1  refer  to  initial  conditions  and  conditions  at  end  of  altitude  interval  dz. 


In -actually  applying  the  above  equations  Fujiwhara  was  forced  to  make  several  further  simplifications 
and  assumptions. 

1.  The  x-axis  waschosen  along  a  number  of  bearings  in  succession  so  that  the  wind  effect  v  sin  ^ 
was  negligible.  This  determined  a  maximum  and  minimum  x  and  defined  the  limits  of  audibility  (See  Figure 
3). 

2.  For  the  observed  atmospheric  structure  an  artificial  one  was  substituted  made  up  of  finite  altitude 
Intervals  in  which  the  temperature  and  wind  gradients  were  linear  (see  Figure  4). 


SIMPLIFIED,  ATMOSPHERIC  STRUCTURE  TO  BE  USED1 
OMPUTiNG!  SOUND*  RAY  PATHS  BY  FUJI  WHAR  A'S  METHOD 
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3.  A -separate  calculation  was  carried  out  for  the  translational  movement  of  the  entire  air  mass 
(consequently  the  sound  ray  as  weE)  for  the.  period  during  which  the  sound  traversed  the  air  mass.* 

”  Tn  1921  E.  A.  Milne  (3)  published  a  theoretical  treatise  on  the  propagation  cf.sound  through  the  atmos¬ 
phere.  By  considering  the  equations  ofrmotionof  a  sound  ray, 


n.  3 


dx  ,  dy  dz  ^ 

gj-  =  lc  +  u,  gi  ---me  +  v,  gp=  nc  +  w 


where:  1,  m  and  mare  direction  cosines  of  normal  towave  surface. at  (x,y,z);  u,  v  ,  w.  are  wind  components 
in  x,  y.and  z  directions  respectively;  c  is  the  velocity  of  sound  due  to. temperature  effects  alone.  Milne  was 
able  to  derive  the  general  solution  for  the  path  of  a  sound  ray  through. the  atmosphere  (see  Appendix  B). 
This  solution  is  given  asv 


IL  4 


1/dl  Sc  ,  du  sv  dw  \ 

1  \dt  qx  ox  B’x  dx  J 

_  1  I  dm  Sc  ©u  »v  0W,\ 

~  4  ay  r  1  6>y  ay  +  n  ayi 

_  i  /  dn  0c  au  .  ^  av  '  aw.y 

"  nr-  \dt  as  1  dz  dz  az  ) 


In  addition  Milne  was  able  to  snow  thafby  assuming  a  stratified  atmosphere  (i.e.  zero  horizontal  wind 
and  temperature  gradients) 


E.  5 


-  =  constant: 
m 


Equation  H.  5  is  most  important  because  it  states  mathematically  that  in  a  stratified  medium  the  normals  to 
the  wave  front  for  any  particular  sound  ray  remain  parallel  to  a  fixed  plane  which  is  itself  perpendicular  to 
the  planes  of  stratification. 

By  choosing  the  x-axis  (equations  E.  4)  in  such  a  manner  that,  for  a  particular  ray,  the  m  direction 
cosine  is  equal  to  zero,  equations  E.  4  may  be  separated  and  integrated  to  give  the  following  resu.’  t: 


*This  translational  movement  is  separate  from  the  refracting  effect  of  the  wind  since  the  former  is  a  motion 
relative  of  the  fixed  point  of  disturbance  while  the  latter  effect  arises  from  the  presence  of  a  vertical  wind 
gradient. 
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It.  6  c  sec  a  +  u  =  c  sec  e  +  u  =  constant, 

o  o 

where  e  =  angle  between  sound  ray  normal  and  the  x-axis. 

Equation'll.  6  Is  the  mathematical  expression  for  the  "Law  of  Sound  Refraction." 

One  more  Investigation  of  Importance  which  has  appeared  recently,  (13)  details  the  work  of  P.  Ruthwell 
and  associates  In  England  in  attempting  a  simple  integration  of  the  equations  of  motion  of  a  sound  ray  In  the 
atmosphere.  The  assumptions  used  in  the  analysis. were  very  similar  to  those  made  by  Fujiwhara  and  Milne, 
namely: 

1.  The  sound  ray  path  was  constrained  to  be  entirely  within  a  vertical  plane. 

2. -  The  atmosphere  was  considered  to  be  stratified  In  horizontal  layers  (no  vertical  wind  currents) 
and  the  observed  structure  could  be  broken  into  layers  where  the  wind  and  temperature  gradients  were  linear 
(see  Figure  4). 

Considering  a  small  altitude  interval  A  y  (Fig..  4),. the  horizontal  component  distance  x*  traveled  by 
a  sound  ray  in  that  interval  and  the  time  required  to  travel  this  distance  are  given  by. 

n.  7  aye  ..  r^B1.(E1(2=%Ri»1/2] 

‘VR.  l-c + r2  -  (r2  (25 + ryjV?  J 

n.  8  t  .  r^~|n2  (25  +B2»1/2  -  (R,  (25  +E!))1M 

2  1  J 

where  the  quantities  c,  R^,  Rg  and  y  are  defined  by  figure  5  below. 

In  the  case  where  R^  =  Rg  (zero  gradient  through  the  altitude  increment  A  y)  equations  H.  7  and 
IL  8  are  indeterminate  and  are  replaced  by 


n.  71 


_ Aye 

(Rj  (2c  +  R1))V 


E.  81 


t  = 


4y  (R2  +c) 
c(Rj  (2c  +  R1))i7^ 


The  additional  horizontal  component  of  sound  ray  *~awel,  x'  *,  due  to  the  translational  movement  of 
the  air  mass  itself  is  given- by 
n.  9  xM=ut 


where,  u  =  mean  wind  velocity  in  altitude  interval  A  y. 


AY  =  small  altitude  interval  through  whiclr  c',  velocity  of  sound  due  to 
temperature,  and  u,  horizontal  wind  velocity  may  be  considered  to 
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The  derivations- of  equations  Tr.  ?  through-IL  9  are  given  in  Appendix  G. 

It  is  obvious  that  a  summation  of  x',  x-  *  and  I  over  a  series  of  altitude  intervals  ( 4  y’s)  will  describe 
the  path  of  a  sound  ray  through  the  atmosphere. 

As  stated  in-the  introduction  to  this  paper  the  geographical-limits  of  the  test  range-for  the  firing  of 
guns  and  explosive  charges  were  such  as  tc  insure  no  blast  damage  to  objects  without  the  confines  of  this 
area.  This  statement  is  true  even  if  the  refraction  effects  of  the  atmosphere  are  considered.  I  bus  if  pro¬ 
pagation  velocity  gradients  in  the  atmosphere  are  such  that  a  blast  ray  (sound  i  ay)  or  a  section  of  a  olast 
wave  is  refracted. downwards  sc  as  to  strike  the  earth  at  a  point  distant  from  ihe  source  of  disturbance,  the 
normal  decrease  in  intensity  (inverse  square  law)  with  distance  assures  that  no  damaging  results  will  occur. 
However,  it  has  been  noted  on  many  occasions  that  damage  has  occured  at  ah  .unreasonably  great  distance 
from  the  explosive  source  and  it.  is  this  latter  phenomenon  which  must  be  investigated. 

Consider  first  the  atmospheric  structure  shown  in  Figure  6a.  For  the  sake  of  simplicity  the  propa¬ 
gation  velocity  illustrated  increases  linearly  with  increase  in  altitude.  From  the  law  of  sound  refraction.il 
is  obviousthat  rays  of  sound  will  be  refracted  downwards  in  such  an  atmosphere.  Figure  6b  shows.the  path 
of  travel  for  four  rays  in  this  atmosphere.  It  will'be  noted  that  rays  1,  2,  3  are  returned  to  different  points 
along  the  abscissa.  However,  because  of  the  change  in  velocity,  gradient  above  the  path  of  ray  3,  ray  4  is 
refracted -in  such  a  manner  that  it  returns 'to  the  same-point  on  the  abscissa  as  ray  3.  Now  if.  the  time  of 
travel  for  rays  3  and  4  in  arriving  at  the  common  point  on  the  abscissa  is  the  same,  constructive- inter¬ 
ference  wiU  occur  and  an  .irregular  pressure  jump  will  take  place.  This  phenomenon  shall  hereafter  be 
termed  "destructive  focusing". 

Returning  now  to  Rothwell’s  equations  (IL  7  to  IL  9)  and  defining  (x*  +  x*  (x*  +  x*  ').,  t„  and  t . 
as  the  horizontal  distance  and  time  of  travel  corresponding  to  the  maximum  altitude  attained  by  rays  3  and 
4  (Fig  6b)  respectively,,  it  is  obvious  that  destructive  focusing  will  occur -when 

IL  10  (x*  +  x'  *)3  =  (x*  +  x*  *)4 


Because  of  the  normal  complexity  of  the  atmospheric  structure  and  the  assumptions  made  in  obtaining 
the  integrated  equations  of  motion,  it  would  be  extremely  difficult  to  obtain  theoretical  solutions  for  equations 
H.  10.  It  is  much  simpler,  and  sufficiently  valid  if  enough  samples  are  taken,  to  obtain  an  approximate  solu¬ 
tion  by  considering  various  simple  atmospheric  structures  and  determining  what  conditions  are  necessary  if 
destructive  focusing  is  to  take  place. 

INVESTIGATION 

As  pointed  out  in  the  preceeding  section,  the  complexities  normally  encountered  in  atmospheric 
structures  prevents  a  simple  solution  for  equations  H.  10.  The  most  obvious  line  of  investigation  points 
tc  numerical  examples  and  subsequent  determination  of  the  factors  necessary  to  produce  destructive 
focusing. 
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For  the  investigation  three  different  atmospheric  structures  were  studied  (see  Figure  7):  positive, 
negative  and  zero  gradients  in  temperature  as  a  function  of  altitude.  The  horizontal  distances  of  travel  and 
the  time  of  travel  were  calculated  fob  several  different  rays  of  sound  in  each  atmosphere.  The  4  y  intervals 
(Figure  7)  were  chosen  so  that  a  particular  ray  of  sound,  ip^,  would  become  horizontal  (i.e.  totally  refracted) 
at  an  altitude  equal  to  A  y.  The  atmospheric  conditions  in  the  layer  above  A  y  necessary  to  focus  destructive¬ 
ly  a  second  ray,  with  tp^  were  then  determined  by  Rothwell-s  equations  and  equations  fL  10. 

As  an  example  consider  the  atmospheric  structure  shown  in  figure  ?a.  The  atmosphere  is  isothermal 
with  a  temperature  corresponding  to  a  sound  velocity  of  HOOTt./sec.  at  the  ground  surface;  the  wind  com¬ 
ponent  along  the  direction  of  travel  is  -10  ft/sec.  If  the  propagation  velocity  increases  with  altitude  as  is 
shown,  it  is  found  that.the  sound  ray,  the  initial  angle  with  the  horizontal  (s^oj  cf  which  is  2°,  will  become 
horizo.ntal  at  an  altitude  of  165  ft.  The  problem  then  is  to  consider  the  wind  and  temperature  structure  of 
the  atmosphere  which  will -focus  destructively  (with  e^o)  sound  rays  e^o,  9^0,.©^*  and  e^o  (i.e.  rays,  the  initial 
anglesof  which  were  4°,  6°,  8°  and  10°  respectively).  Thus  if  x^  t^  and  x,j,  tg  represent  the  horizontal  dis¬ 
tance  and  time  of  travel  of  rays  o9°  and  9^0  through  the  interval  A  y,  (165  ft.  in  this:case)  and  if  x^,  t^  repre¬ 
sent  the  same  quantities  for  ray  e^  in  altitude  interval  A  y^  above  a  y^,  the  conditions  necessary  for  destruc¬ 
tive  focusing  of  rays  «9<>  and  0^°  are 
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xv,  x0,  t  and  t„  are  easily  solved  from -the  numerical  .values  given  after  equations  EL  2.  Since  JR,  Is  aiso 

1  u  1  '6  D 

known,  equation  EL  3  may  be  solved  tc  flnd  the  value  of  Cg.  Once  this  Is  done  substitution1  of  c^  Into  either 
equation  EL.23  or  equation  IIL  26  will  result  in  a  solution  for  Ay^  and  the  atmospheric  structure  will  be 
determined.  In  the  numerical  example  above  A  is  found  to  be  209.0  ft.  and  equal  to  1100.  15  ft./sec. 

The  solution  for  Cgin  equatlonEL  3  may.be  accomplished,  without  too  much  difficulty,  by  graphical 
methods.  However,  another  solution  may  .be  found  (see  Appendix  D)  for  the  quantities  A  y2  and  by  the 
following  formulae:  -  - 


HI.  41 


;  ,  R5  , 

2  cos  hHL  _i  „ 
**2 


m  » 


2(R5)  (X1  "  X2^ 


The  equations  IIL  41  and  EL  42,  while  exact,  necessitate  a  determination  of  t^  and  tg  accurate  to  five 
significant  figures.  Unfortunately  the  simplified  equations  of  RothweE’s  are  not  always  accurate  to  I’*is 
number  of  significant  figures  so  that  results  obtalned.in  solving  for  Cg  and  a  y2  by  means  of  equations  EL  41 
and  EL  42  may- oftentimes  be  totally  Invalid.  For  this  reason  the  graphical  solution  for  equation  EL  3- was 
adhered  to  In  aH  calculations. 

Details  of  the  calculations  Involved  in  studying  the  atmospheric  structures  of  Figure- 7  wer  e  not  con¬ 
sidered  essential  for  the  main  body  of  the  report  and  consequently  these  may  be  found  at  the  end  of  the  re- 
port-ln  Appendix  E. 

CONCLUSIONS 

The  last  section  dealt  very  briefly  with  the  method  used  In  Investigating  those  atmospheric  conditions 
wliich  are  conducive  to  destructive  focusing.  The  description  was  purposely  kept  brief  because  of  the  relative 
unimportance  of  the  mathematical  gymnastics.  The  numerical  results  for  one  computation  were  considered 
sufficient  to  IEustrate  the  method. 

The  systematic  Investigations  of  the  three  ideal  atmospheres  pictured  In  Figure  7  revealed  one  very 
significant  and  consistent  result.  In  order  that  destructive  focusing  of  two  or  more  sound  rays  might  occur, 
a  temperature  inversion  in  the  atmosphere  was  found  necessary.  Referring  to  Figures  7a  and  7c,  a  positive 
temperature  gradient  above  the  altitude  Increment  Ay^  was  a  necessary  condition  for  destructive  focusing 
to  take  place.  For  the  atmosphere  illustrated  in  Figure  7b,  a  negative  or  positive  (or  zero)  temperature 
gradient  above  the  altitude  Interval  A  was  a  necessary  condition.  Actually  the  results  from  the  study  of 
Figure  7b  might-  sound  ambiguous,  but  by  definition  a  temperature  inversion  in  the  atmosphere  Indicates 
merely  the  existence  of  a  positive  temperature  gradient.  K  one  particular  line  of  propagation  (direction 
from  the  source)  is  considered  —  such  as  the  examples  in  Figure  7  —  certain  conditions  must  also  be  fulfilled 
in  the  wind  gradient  structure.  However,  the  problem  is  concerned  with  destructive  blast  propagation  in  any 
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direction  and  consequently  the  wind  component,  Uj  may  vary  infinitely  for  each  atmospheric  structure  existent 
Ih  nature.  Thus  it  is  only  practicable  to  determine  the  necessary  rather  than  the  sufficient  condition. 

Earlier  in  the  paper,  mention  was  made  of  the  fact  that  water  vapor  present, in  the  atmosphere  also 
has  an  effect  on  sound  propagation.  Since  this  element  (water  vapor)  varies  only  with  altitude  —  neglecting 
local  and  low  lever  influences,  such  as  a  body  of  water  —  its  effect  may  be  considered  simultaneously  with 
the  temperature  gradient.  In  Appendix  A,  formulae  are  derived  for  the  propagation  of  sound  in  dry  and  moist 
air.  Figure*8  below  illustrates  graphically  the  effect  of  water  vapor  on  the  velocity  of  sound  at  sea  level  and 
and  at  3000  feet  if  altitude.  Thus  In  constructing  the  curve,  showing  the  variation  of  sound  velocity  with  alti¬ 
tude  as  a  function  of.  temperature  alone  it  is  necessary  to  consider  the  correction  due  to  the  presence  of 
water  vapor  before  the  true  picture  can  be  given  of  variation  of  propagation  velocity  with  altitude.  As  can 
be  seen  from  figure  8  the  mere  presence  of  water  vapor  in  the  atmosphere  would  not  seriously  affect  the 
temperature-propagation  curve,  but  a  vertical  gradient  in  the  relative  humidity  might  possibly  change  the 
temperature-propagation  curve  from  one  of  a  negative  or  zero  gradient  to  a  positive  gradient  and  thus  bring 
about  a  sufficient  condition  for  destructive  focusing. 

Although  the  existence  of  a  positive  temperature  gradient  may  be.  absolutely  determined  ohlyby  an 
actual  sounding  of  the  atmosphere,  there  are  many  meteorological  phenomena,  observable  from  the  ground, 
which  indicate  its  presence.  The  mere  fact  that  the  normal  or.  average  temperature  structure  in  the  atmos¬ 
phere  is  characterized  by  a  negative  gradient  strongly  suggests  the  association  of  abnormal  gradients  with 
observablembnormalties..  A  discussion  of  pertinent  meteorological  elements  is  given  below. 

T.  Wind  effects.  A  strong  horizontal  wind  tends  to  produce  considerable  eddy  currents. either  near 
large  .physical  objects  on  ground  surfaces  or' along  shear,  surfaces  in  the  atmosphere  itself.  These  eddy 
effects  thoroughly  mix  the  atmospheric  elements  and  produce  the  normal- negative  temperature  gradient. 

The  same  eddy  effects  produce  vertical  wind  currents  and  would  distort  sound  wave  fronts.in  an  unpre¬ 
dictable  inanner.  Since  the  entire  theory  of  refraction  in  the  atmosphere  of  sound  or  blast  waves  is  based 
upon  the  assumption  of  laminar  flow  with  no  vertical  currents,  the  eddy  effects  produced  by  the  strong  hori¬ 
zontal  flow  of  air  would  nullify  any  results  of. the  above  investigation.  However,  the  distortion  of  sound 
waves  produced  by  such  a  flow  would  make  destructive  fpcusing  extremely  improbable  and  the  presence  of 
strong  winds  reasonably  indicates  conditions  not  conducive  to  destructive  focusing.  Strong  surface  winds 
aloft,  and,  although  conclusive  evidence  is  not  at  the  moment  available,  surface  winds  of  15  mph  or  greater 
seem  likely  meteorologically  safe  firing  conditions. 

2,  Temperature  Inversions  A  temperature  inversion  in  the  atmosphere  accompanies  numerous 
synoptic  weather  phenomena  but  most  of  these  phenomena  aie  associated  with  a  rather  special  synoptic 
pattern.  In  the  preceding  paragraph  the  assertion  was  made  that  a  strong  horizontal  air  flow  tended  to 
produce  a,  negative  temperature  gradient  in  the  atmosphere.  Consequently,  if  temperature  inversions 
are  to  exist  the  horizontal  transport- of  air  must  be  small.  In  general,  light  winds  —  up  to  any  appreciable 
altitude  —  are  found  within  high  pressure  areas  which  tend  to  stagnate  for  periods  of  several  da  s’  durat¬ 
ion  over  a  large  ground  area.  If  the  high  pressure  cell  on  the  earth’s  surface  is  small,  say  of  the  order 
of  100  miles  in  diameter,  and  is  moving  fairly  rapidly,  the  winds  although  light  or  calm  at  the  surface  will 
be  found  to  be  high  enough  at  2000-3000  feet  in  altitude  to  produce  mixing  and  consequently  a  negative 
temperature  gradient  at  that  altitude.  True,  close  to  the  earth's  surface,  there  may  be  a  temperature 
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Inversion  due  to  the  greater  radiation  of  the  ground  (compared  to  the  air)  but  this  Inversion  would  be 
sufficiently  low  to  keep  destructive  effects  wiihln  the  blast  range  area.  The  large  scale  high  pressure  area 
which  stagnates  for  several  days  over  portions  of  the  earth’s  surface  is  characterized  by  the  -ubsld’.ng  .r 
sinKing  of  the  air  mass  and  ^corresponding  increaoe  in  the  temperature  gradient  (le~s  negative  if  original¬ 
ly  negative  and  more  positive  if  ot  iglnally  positive),  caused  by  the  adiabatic  heating  ir.  compression.  The 
wind. velocities  are  fairly  low  up  to  altitudes  of  perhaps  600C  or  7000  feet  and  the  atmosphere  becomes  quite 
stratified.  This  condition  is  quite  conducive  to  destructive  focusing  and  there  is  little  distortion  of,  sound 
wave  fronts  in  passing  through  t'ne  atmosphere.  These  extensive  areas  ~f  high  pressure  tend  tc  be. more 
frequent  and  more  perslsten;  in  the  summer  and  wither  months  and  they  may  easily  be  recognized  either 
by  reference  to  the  current  weather  map  or  by  observing  meteorological  elements  such  as  temperatures, 
visibilities  and  cloud  types. 

There  exist  weather  situations  in  which  temperature  inversions  occur  in  low  pressure  areas  in  the 
presence  of  fairly  high  wind  velocities.  The  presence  of  warm  or  c~ld  frontal  surfaces  in  the  immediate 
vicinity  indicate  temperature  inversions  at  the  boundry  between  the  different  air  masses.  However,  the 
slope  of  the  frontal  surface  implies  a  horizontal  temperature  gradient  as  well  as  a  vertical  one  and  the 
laws  of  sound  refraction  again  do  not  apply.  Of  course  the  immediate  presence  of  a  -warm  or  cold  front 
usualiy  indicates  rather  adverse  weather  conditions  —  rain,  snow,  high  winds,  poor  visibility,  etc.  —  and 
the  relatively  rapid  movement  of  the  frontal.surface,  in  a  horizontal  direction,  would  make  destructive 
focusing  rather  difficult  in  view  of  the  wind  discontinuity  at  the  frontal  surface. 

Ultimately  the  decision  as  to  whether  or  not  a  firing  is  to  take  place  should  be  based  upon  a 
sounding  of  the  atmosphere  by  suitable  meteorological  instruments.  However,  a  perusal  of  the  preceeding 
paragraphs  enables  one  to  make  some  good  guesses  as  to  the  possibility  cf  destructive  focusing  for  a 
given  atmospheric  situation.  Some  rules  of  thumb  might  be  expressed  as  follows: 

1.  A  high  pressure  area  which  is  stationary  over  the  area  in  question  becomes  dangerous  after 
the  first  day  of  its  presence. 

2.  Large  diurnal  temperature  variations  at  the  earth’s  surface  indicate  unfavorable  firing  conditions. 

3.  Poor  visibilities  and  light  winds  in  the  early  morning  followed  by  cloudless  conditions  throughout 
the  day  indicate  dangerous  focusing  possibilities. 

4.  Poor  visibilities  and  light  winds  in  the  early  morning  followed  by  the  formation  of  cumulus 
(swelling)  clouds  in  the  late  morning  and  early  afternoon  indicate  safe  conditions  for  an  afternoon  firing. 

5.  Light  surface  winds  and  low  relative  humidities  at  the  ground  plus  the  presence  of  a  stratified 
cloud  deck  below  10,000  feet  indicate  dangerous  conditions  for  firing. 

When  soundings  of-the  atmosphere  are  to  be  taken  a  rather  important  question  arises  as  to  how  high 
or  to  what  altitude  the  atmosphere  must  be  examined.  This  question  cannot  be  answered  simply  because  cf 
the  dependence  of  the  limits  of  destruction  upon  the  weight  of  the  detonating  charge  (See  Sachs’  L  1).  How¬ 
ever,  by  means  of  equation  IL  7  one  may  plot  x*  as  a  function  of  Ay  and  (Rg  -  R^)  and  assuming  various 
values  for  the  quantity  c.  It  is  foun  that  over  a  normal  range  of  values,  variations  in  c  have  little  effect 
on  the  dependent  variable  x’  and  may  be  considered  constant  for  our  purposes.  Now  if  Rg  is  set  equal  to 
0,  x*  will  represent  the  horizontal  distance  traversed  by  a  sound  ray  from  the  origin  (source  of  disturbance) 
to  the  point  where  the  noimal  to  the  wave  front  becomes  horizontal  and  thus  2x’  would  define  the  point  at 
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which  the  sound  ray  returns  to  the  earth's  surface  (neglecting  earth  curvature  effects).  Figure  9  below  re¬ 
presents  the  plot  of  equation  EL  7.  The  physical  significance  of  such  a  graph  Is  that  if  the  actual  velocity 
gradient  of  sound  existent  in  the  atmosphere  is  superimposed  upon  the  graph,  only  that  portion  of  the 
gradient  curve  wliich  lies  between  the  abscissa  and  the  parametric  cuive  2x*  =  n  need  be  considered  in 
determining  whether  or  not  destructive  focusing  will  occur  between  the  origin  and  and  2x*  =  n.  For  example, 
supposing  the  propagation. velocity  gradient  along  a  certain  direction  from  a  detonating  source  is  represented 
by  curve  in  Figure  9.  For  a  particular  weight  of  explosive  we  find  by  Sach's  formula,  L  1,  that  the  destruc¬ 
tive  effects  will  be  confined  to  a  circular  area  of,  say,  10  mile  radius  (i.e.  2x*  =  10  miles).  In  Figure  9,  the 
gradient  curve  crosses  the  parametric  curve  2x*  =  10  at  an  altitude  of  3000  feet.  Therefore,  only  the  region 
of  the  atmosphere  below  3000 ‘feet  need  be  considered  for  determining  the  existence  of  destructive  focusing. 
By  means  of  this  plot  and  the  meteorological  data  on  vertical  velocity  gradients  existing  12  hours  or  less 
before  the  planned  firing  time  (available  from  local  weather  stations),  a  rough. estimate  of  the  height  to  which 
the  atmosphere  must  be  sounded  (to  determine  the  existence  of  destructive  focusing)  is  obtained. 

The  preceedlng  paragraphs  have  been  concerned,  mainly  with  implications  rather  than  objective  facts 
pertaining  to  the  problem  of  destructive  focusing.  It  is  felt  that  such  a  superficial  discussion  is  permissible 
within  the  scope  intendedior  such  a  preliminary  investigation.  However,  an  attempt  was  made  to  prove  some 
of  the  conclusions  stated  in  earlier  paragraphs  by  examining  meteorological  conditions  on  days  when  fairly 
valid  complaints  were  received  from  inhabitants  living  outside  the  Proving  Ground  area.  Three  examples 
of  typical  complaints  were  chosen  at  random  (February  20,  1944,  March  16,  1944  and  March  26,  1945)  and 
the  atmospheric,  conditions  examined.  In  each  case- the  following  conditions  were  noted: 

1.  A  positive  temperature  gradient  existed  at  some  point' in  the  atmosphere  (at  relatively  low 
altitudes).. 

2.  Surface  winds  were  less  than  8  mph. 

3.  Visibilities  were  poor  in  either  fog  or  smoke. 

4.  A  high  pressure  area  (surface)  characterized  the  synoptic  situation. 

This  evidence,  while  far  from  conclusive,  does  certainly  point  to  atmospheric  conditions  conducive  to 
destructive  focusing. 

On  January  14,  1948,  a  very  clear  cut  case  of  destructive  focusing  was  encountered  and  accurately 
reported.  In  this  instance  a  12,000  lb.  charge  was  detonated  in  the  face  of  a  surface  wind  blowing  from  the 
north  at  approximately  14  mph.  At  4500  yards  to  the  north  only  a  faint  sound  was  heard;  15,000  yards  to  the 
north  no  sound  was  heard;  25,000  yards  to  the  north  a  sharp  crack  followed  by  a  prolonged  rumbling  was 
noted  by  numerous  observers.  The  normal  decrease  in  intensity  of  the  blast  wave  moving  up  wind  was 
quite,  in  line  v/ith  the  observations  at  4500  and  15,000  yards  but  some  change  in  atmospheric  structure 
aloft  was  necessary  to  produce  the  observed  effects  at  25,000  yards.  Figure  10  illustrates  the  sound  pro¬ 
pagation  gradient  along  the  line  from  the  detonating  source  to  the  point  at  25,000  yards  where  the  report  was 
heard.  From  the  figure  it  is  eas’ly  seen  that  the  necessary  condition  for  destructive  focusing  existed. 

One  final  comment  should  be  added  at  this  point.  It  has  been  suggested  that  the  destructive  effects 
of  large  detonating  charges  might  be  attributed  to  shock  transmitted  by  the  ground.  There  are  two  facts 
which  strongly  discredit  this  possibility: 

1.  If  the  transmitted  shock  wire  a  surface  v/ave  (in  the  earth's  crust)  it  would  decrease  approxi- 
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mately  as  the  Inverse  square  distance  from  the  source  and  residents  living  close  to  the  Proving  Ground 
.would  suffer  far  in  excess  of  those  living  at  some  distance-from  the  area  —  a  fact  not-observed  In  actual 
cases. 

2.  K  the  shock  wave  were  reflected  or  refracted  from  some  Interlace  or  change  in  structure  below 
the  earth's  surface,  a  particular  location  would  suffer  equally  from  every  charge  (of  equal  magnitude) 
detonated  at  the- same  point  on  the  Proving  Ground  —  likewise  a  fact  not  observed  during  firings; 


W.  Berning 
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APPENDIX  A 

Derivation  of.tha  formulae  for  the  propagation  of  sound 
in  the  atmosphere  (References  (14)  and  (15) ) 

Consider  an  . Imaginary  tube  of  unlt  cross  section  parallel  to  the  line  of  propagation  and  containing 
a  portion  of  the  sound  wave  front.  Let  -F-  be  the  external  force  supplied  so  as  to  keep  the  form  of  the  wave 
front  unchanged  and  moving  with  a  uniform  velocity  a.  The  momentum  entering  one  end  of  the  imaginary 
tube  per  interval  of  time  is  equal  to  the  momentum  leaving' the  tube  during  the  same  time  Interval.  Thus 
if  P  represents  the  external  atmospheric  pressure,  AP  the  pressure  Increase  in  the  sound  wave,  a,  the 
velocity  of  entry,  of  air  of  density  pand  a-+  u  and  p^  the  velocity  and  density  of  the  air  as  it  leaves  the  tube 


(A.  1) 

M  (a  *+  u)  +  P  +  A  P 

or  _ 

(A.  2) 

p  au=P  +  F. 

Now,  if  F  is  equal  to  0,  the  case  which  corresponds  very,  closely  to  sound  propagation  in  the  atmos¬ 
phere  since  the  pressure  behind  the  sound  wave  is  almost  exactly  equal  to  that  ahead  of  it 


(A.  3) 

pa  =  AP 

u 

or 

(A.  4) 

a=  AP 
Pu 

However,  from  the  statement  of  the  problem 


(A.  5) 


Pj  (a  +  u)  =  Pa 


and  if  the  ccmpressional  changes  take  place  adtabatically  p^  is  given  by 


(A.  6) 


P+  AP 
P  ' 


y  -  ratio  of  specific  heats  in  air 


so  that  equation  (A.  4)  may  be  written 


(A.  7) 


if  we  neglect  the  product  A  P.  u. 

Equation  (A.  7)  is  the  general  expression  for  the  velocity  of  sound  in  a  dry  atmosphere  where  the 
density  P  of  the  atmosphere  is  dependent  only  upon  temperature  and  pressure. 
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In  practically  treating  the  problem  of  sound  propagation  through  the  atmosphere,  account  must  be 
taken  of  the  presence  of  water  vapor.  Because  of  Its  composition  water  vapor  has  a  lower  density  than  the 
air  which  It  replaces  ahd  consequently  moist  air  Is  lighter  than  dry  air.  The  equation  of  state  for  a  perfect 
gas  Is  given  by 

T  »  absolute  temp. 

P  -  density 
P  «  pressure 
M  =  molecular  wt.  of  gas 
r*  =  universal  gas  const. 

Dalton’s  Law  states  that  the  sum  of  the  partial  pressures  of  a  mixture  of  gases  Is  equal  to  the  total  pres¬ 
sure.  Thus  If  Pm  Is  the  density  of  a  mixture  of  gases  of  densities  P  ^  (dry  air)  and  p  w  (water  vapor),  the 
equation  of  state  (A.  8)  becomes 


(A.  9) 


D  -  M  (p-g)  .  M  w  ( e) 
m  ~  R*T  +*R*T 


M  w  =  molecular  wt.  of  water  vapor 


PM/,  0.379  <  \  €  =  partial  pressure  of  water  jpor 

~R*T  \  PI 


Substitution  of  p  into  equation  A.  7  and  designating  a  and  a*  as  the.veloclty  of  sound-in  dry  and  moist  air 
respectively,  one  obtains 


(A.  10) 


If  now,  the  medium  In  which  the  sound  wave  is  being  propagated  is  Itself  moving  relative  to  a 
coordinate  reference  system  (such  as  the  earth)  the  absolute  velocity  (with  respect  to  the  earth)  of  the 
sound  wave,  designated  as  a"  is  given  simply  by 


(A.  11) 


R=  wind  velocity  normal  to  wave  front, 


& 

c 

y> 

f  . 


f 
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APPENDIX  B 

Condensed  derivation  of  the  "Law  of  Sound  Refraction"  as  taken  from  Milne’s,  work  (Ref.  (3) ) 


Let:  u,  v,  w  be  the  components  of  velocity  at  point  (x,  y,  z)  in  the  atmosphere  moving  with  a  con¬ 
tinuous  flow;  a  be  the  velocity  of  sound,  a  function  of  x,  y,  z;  1,  m,  n  De  the  direction  cosines  of  the  normal 
to  the  wave  front  at  point  (x,  y,  z).  Then  the  equations  of  a  sound  ray  are  given  by 


B.  1 


dx 

ar 


=  la  +  u, 


dy 

dt 


=  m  a  +  v. 


^!  =  na  +  w 
dt 


where  t  is  the ’time. 

The  ray,  however,  is  not  simply  obtained  by  a  simple  integration  of  equations  B.,  1  since  1,  m  and  n 
are  known  only  at  their  initial  values.  Therefore  let  (x,  y,  z)  and  (x  +j  x,  y  +  5  y,  z +  .5  z)  represent  two 
points  of  intersection  of  neighboring  sound  rays  with  one  sound  wave  front  at  time  t.  Subtracting  pairs  of 
equations  of  type  B.  1  one  obtains, 


=  6  (la  +  u)=£dxJ~-(la  +  u) 

B.  2  =  6  (ma  t  v)  y  -5-y-(ma  +  v) 

---J;-  =  5  (na  +  w)  =£  dx~fr^na  +  w) 


Now  differentiating  the  relation 


B.  3 


1  fix  +  m  6  y  +  n{z  =  0 


I 

1 

4 

f 


with  respect  to  time  and  substituting  from  B.  2  for  d  {  d  x)/dt,  etc.,  one  obtains,  after  reduction, 


B.  4 


m 

ax 


+  m 


av 

ax 


+  n 


aw  \ 
ax  j 
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1  - 
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Equation  3.  4  holds  for  all  values  of  {  x;  d  y;  6  z  satisfying  B.  3,  hence. 


B.  5 


r( 


dl 

dt 


5a  6u_ 
9x+1  gx 


+■  m 


dv 

d* 


-f  n 


dw\ 

ax j 


/  dm 

\  3T  + 

(«*  + 

\dt 


^  +1 

ay 


6a 

ez 


+  i 


au  av 

ay  +m  ay  +n 


au 

ez 


+  m 


av 

ez 


+  n 


aw\ 

ay/ 

aw\ 
az ) 


Equations  B.  1  and  B.  5  give  the  fay  completely.  Equation  S.  5  may  be  expressed  more  simply  as  follows. 
Let 


V  =  a  +  lu+mv  +  nw 


and  let  l,*m  and  n  be  considered  as  constant  in  differentiating  with  regard  to  x,  y,  z.  Then  B.  5  may  be 
written 


1 

f 


j 


* 


i 


B.  6 


1  /dl 
1  Vdt 


=  V 


)- 

(* 


1  /dm  ev\_  1  /  dn 
m  \  "3r  ay./  n  '  3F+ 


ex 


+  m 


av\ 

ez/ 


or  explicitly, 


B.  7 


dl 

dt 


(2  2  ’j 

^m  +n  ) 


0V 

ax 


.  av  . 

+ 1  m  —  +  1-n 

ay 


evv 
az  j 


Equations  B.  1  express  the  direct  propagation  of  the  sound,  as  influenced  by  the  bodily  displacement  of  the 
medium,  while  equations  B.  6  express  the  refraction  of  the  sound  caused  by  the  variation  of  the  normal 
speed  V  over  the  surface  of  the  wave  front. 
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If,  as  la  this  paper,  we  are  concerned  only  with  the  motion  taking  place  within  a  stratified  medium,, 
the  z-axlscan  always  be  chosen  so  that  the  velocity  component  w  =  0  and  u,  v,  a  are  functions  of  z  only.  In. 
this  case  equations  B.  1  and  B.  6  reduce  to 

„  n  dx  ,  dy  dz 

B;  8  gj-  =  la  +  u,  j  =  ma  +  V|  gj-  -  n  a 


B.  9 


1  dl  1 
rar  =  m 


dm' 

a r! 


/dn 

lar+ 


=  n 


av 

az 


V  =  a  + 1  u  +  m  v 


The  first  of  equations  B.  9  Integrates  immediately  In  the  form 


Bi.lO 


—  =  constant, 
m 


Equation  B..10  states  that  in  a  stratified  medium  the  normals  to  the  wave  fronts  along  a  particular  sound 
ray  remain  parallel  to  a  fixed  plane  which  is  itself  perpendicular  to  the  planes  of  stratification  (in  this 
instance  the  planes  parallel  to  the  xy  plane). 

If  now  one  chooses  the  x,  y  axes  so  that  the  projections  of  the  normals  are  parallel  to  the  x  axis, 
m  =  0  and  so 


dz 

dt 


=  n  a, 


i-M 


aa 

62 


+  1 


6u  \ 

dz)  ' 


whence, 


—  --L  /  £1+1  £1  \ 

dz  a  \  dz  6z  / 


or, 


B.  11 


d  / a\  du 

"•aruraa 
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APPENDIX  C 


Derivation  of  Rothwell's  equations  (see  reference  (13))  for  the  motion  of  a  sound  ray  in  the  atmosphere. 


The  "refraction  law'1  states  that 


c  sec  e  +  u  =  A, 


e  =  velocity  of  sound  in  calm  air 
u  =  horizontal  wind  velocity  component 
A  =  constant 

e  —angle  between  horizontal  and  normal  to  the 
wave  front. 


This  equation  may  be  written  in  the  form 
sec  ©  -  (R/c)  +  1, 


where 

C.  3  R  =  A  -  (c  +  u). 

Figure  C— 1  represents  a  "linearized"  vertical  segment  of  the  atmosphere  (explained  in  the  text).  In 
Figure  6—1  let  u,  c  and  R  at  the  layer  boundaries  be  distinguished  by  subscripts  1  and  2,  and  let 
C.  4  c  -  c  =  hy 

C.  5  Uj  -  u  =  ky 

y  being- measured  from  boundary  1.  Ay  Is  the  thickness  of  the  layer. 

From  equations  C.  1,-C.  4  and  C.  5  are  derived  the  following: 


Cj  sec  0  -  (A  -  Uj) 
y  =  h  sec  0  +  k 

dy  + 

as  *  '  - - — j! — 

(h  +  k  cos  e) 

Ak  +  (Ujh  -  c^.  k)  sec  © 
u  h  sec  e  +  k 


C.  3  r*  -  °1  k  +  (A  “  u|)  h 

'  h  sec  e  +  k 

i 

hi  Figure  C-2  which  represents  a  ray  PQ’  in  an  elementary  layer  of  infinitesimal  thickness  gy,  PQ 
is  the  direction  of  the  wave  normal  and  QQ‘  is  the  component  wind  displacement  of  an  element  of  wave  front 
while  traversing  the  layer.  The  velocity  in  the  layer  is  c,  hence  the  time  6 1  taken  by  the-  element  to  cross 
the  layer  is 

6  t~  g  y/c  sin  e. 


6  x  =  d  x*  -  6  x"  =  g  y  cot  e  +  u  6 1. 
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APPENDIX  D 


to  Figure  D-l,  Ay  represents  the  altitude  at  whlchthe  sound  ray.(l)  becomes  horizontal.  The 
problem  is  to  determine  the  mean  temperature  c2  and  the  altitude  Ay2,  above  A  y^,  for  a  second  ray  (2). 
which  will  destructively  focus  with  ray  (1).  Ifx^,  t.,  and  Xg.  tg,. represent,  respectively,  the  horizontal 
distance  and  time  of  travel  of  sound  rays  (1)  and  (2)  in  the  altitude  interval  A  y^,  and  if  x^,  t^  represent 
the  horizontal  distance  and  time  of  travel  of  sound  ray  (3)  in  altitude  interval  ^y,,,  the  conditions  for 
destructive  focusing  (at  the. ground)  of- these  two  rays  and  expressed  by 


AH  t3  =  V-t2 


D.  12  x3  =  xi  “  x2 


The  equations  to  be  solved  then  are  (See  IL  7,  IL  8) 


D  ,  21  t  =  ^ 

-52E; 


D.  22  2A  _R1 

3  2 


,/2} 

K<:!c„*r.))  1/2t 


where  A,  R^,  and  R2  (=0)  are  determined  for  a  particular  ray  in  a  particular  atmosphere  (See  Append jc  C). 


From  equation  D.  21 


|RX  (2c2  +  R1)J1/2  13T1 


-  2  2 
C  2  ^2 


2c2  +  Ri  tgRj 


-  4y2  A*2JA72  ■ 

“2=tI^+ 
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whence; 
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or,  after  simplification 


s 
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APPENDIX  £ 

Summary  of  calculations  for  the  three  atmospheric  structures  studied  in  the  text. 

The  equations  used  in  determining  the  path  of  a  sound  ray  through  the  atmosphere  (see  appendix  C) 
are  given  as  - 


c 


E.  3 


t'F]v%)l(^(25tR^V2-(s‘(25+R'))1/2 


Now,  since  the  total  horizontal  range  x  of  a  sound  ray  in  a  segment  of  the  atmosphere,  ^  y,  is  given 
by  x  =  x*  +  x",  equations  E.  1  and  E.  2  may  be  combined  to  give 


E.  4 

X=  Ay 

A  i/o 

in7  S  +  Ri-! 

R^  (2c  +  R^)  ^ 

Rl+R2 

/ 

U  +  R2- 

(B2<2S  +  E2>) 

1 1/2 
1/2 


} 


The  investigation  procedure  for  examining  the  three  idealized  atmospheres  (Figure  7)  was  explained 
adequately  in  the  text  and  only  the  results  will  be  given  in  this  appendix, 

1)  Example  (a)  Figure  7  a. 

For  this  example  an  isothermal  structure  was  considered  and  an  altitude  interval  A  y^,  chosen  so 
that  a  sound  ray  of  initial  angle  e  =  2°  (angle  with  respect  to  horizontal)  would  be  just  refracted  in  the 
interval.  The  horizontal  range  and  time  of  passage  of  the  sound  ray  e  =  2°  in  the  layer  a  y^,  was  then  cal¬ 
culated  by  means  of  equations  E.  3  and  E.  4.  The  horizontal  ranges  and  times  of  passage  of  sound  rays 
e  =  4°,  e  =  6°,  e  =  8°  and  ©  =  10°  were  calculated  for  the  same  altitude  interval  A  y^.  Table  I  gives  the 
values  for  x  and  t  of  sound  rays  e  =  2°,  e  =  4°,  ©  =  6°,  ©  =  8°  and  e  =  10°  in  the  altitude  intervalAy^,  In 
addition  the  various  constants  for  use  in  equations  E,  3  and  E,  4  are  tabulated. 
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TABLET 


Ft.  sec. 

.  Ft/sec. 

Ft. 

Sec. 

Ft. 

Sec; 

Ft. 

0 

ci 

C2 

C 

U1 

“2 

u 

A 

R1 

R2 

Y1 

t 

X 

t  (2) 
-t  (0) 

X(2) 

2° 

1100 

1100 

1100 

-10 

-9.34. 

-9.67 

1090.66 

0.66. 

.0 

165 

8.660 

9445.5 

0 

0 

4° 

iioo 

1100 

1100 

-10 

-9.34 

-9,67 

1092.64 

2.64 

L98 

165 

2.324 

2532.3 

6.336 

6913.2 

6° 

1100 

1100 

1100 

-10  _ 

-9.34 

-9.67 

1096.38 

6.38 

5.72 

165 

1.429 

1548.3 

7.231 

7897.2 

8° 

1100 

1100 

1100 

-10 

-9.34 

-9.67 

1100.78 

10.78 

10.12 

165 

1.096 

1185.7 

7.564 

8259.8 

10° 

1100 

1100 

1100 

-10 

-9.34 

-9.67 

1106.94 

16.94 

16.28 

165 

0;873 

937,7 

7.787 

8507.8 

0) 

(2) 

(3) 

(4) 

(5) 

(6) 

(?) 

(8) 

(9) 

(10) 

(11) 

(12) 

(13) 

(14) 

(15) 

The  equations  for  computing  A,  and  R  are  given  by 

c  sec  0  +  u  =  A, 

A  -  (c  +  u)  =  R. 

In  order  that  destructive  focusing  Is  to  occur  between  0  =  2°  and  0  =  4°,  6  =  6°,  0  =  8°  and  0  =  10°, 

It  Is  obvious  that  the  time  of. travel  and  the  horizontal  range  for  rays  0  =  4°,  6°,  8°  and  10°  In  an  altitude 
Ay2  above  4 y^  must  be  equal  to  the  time  and  distance  differences  (columns .(14)  and  (15)  of  Table  I).between 
0  =  2°  and  the  other  sound  rays.  Substituting  these  time  and  distance  differences  Into  equations  E.  3  and 
E.  4  t  babies  one  to  solve  for  the  atmospherlc  structure  of  altitude  Interval  A  y2-  Table  n  tabulates  the 
characteristics  of  the  atmosphere  above  Ay^  which  will  produce  destructive  focusing. 

TABLE  E 

Ft/sec.  Ft/sec  Ft. 

0  C  Cl  C2  mn1,  u2  Ay2  R1  R2  A 

4°  1100.15  1100  1100.30'  -9.34  -  7.66  209.0  1.98  0  1092.64 

6°'  1106.4  1100  1112.8  -9.34  -15.4  406.2  5.72  0  1096.38 

8°  1109.2  1100  1118.4  -9.34  -17.6  565;4  10.12  0  1100.78 

-  10°  1114.6  1100  1129.2  -9;34  -22.3  739.0  16.28  0  1106.94 

Table  a  Illustrates  quite  clearly  the  temperature  function  necessary  for  destructive  focusing. 

2)  Example  (b)  Figure  7b. 

In  example  (b)  a  slightly  different  procedure  was.  followed.  Three  sound  rays,  0  =  2°,  .0  =  4°,  and 
0  =  6°,  were  traced  through  an  altitude  interval  Ay^  which  was  allowed  to  vary  but  within  which  the  mean 
temperature  was  held  constant.  Effectively  this  means  that  three  separate  atmospheres  were  studied 
(temperature. and  wind  gradients  varied  over  Ay^).  However,  the  wind  gradients  were  so  chosen  that  the 
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sound  ray  0  =  2°  was  just  refracted  In  the  altitude  interval  4  Equations  E.  3  and  E.  4  were  then  used 
to  determine  the  temperature  structure  and  altitude  interval  ^y^  which  would  destructively  focus  sound 
rays  0  =  4°  and  0  =  6°  with  the  ray  0  =  2°.  The  results  of  these  calculations  are  given  in  .Table  m  below. 


6 

Ft/sec 

C1  C2 

C3 

TABLE  HI 

Ft/sec 
51  U1 

*2 

U3 

Ft. 

% 

av2 

52 

2° 

1100 

1110 

1105 

10 

0.66 

— 

200 

0 

— 

4" 

1100 

1110 

1101 

1105 

10 

0.66 

11.64 

200 

253.4 

1105.5 

6“ 

1100 

1110 

1109.9 

1105 

10 

0.66 

6.48 

200 

491.7 

1109.85 

4° 

1100 

1110 

1103.8- 

1105 

10 

0.66 

8.84 

400 

507.2 

1106.9 

6° 

1100 

1110 

1109.9 

1105 

10 

0.66 

6.48 

400 

983.3 

1109.85 

4° 

1100 

1110 

1103.3 

1105 

10 

0.66 

8.84 

600 

760.8 

1106.9 

6° 

HOC 

1110 

1109.9 

1105 

10 

0.66 

6.48 

600 

1474.9 

1109.85 

Cf,  Cg,  and  Ug  represent  the  sound  velocity  and  wind  velocity  at  the  lower  and  upper  limits, 


respectively,  of  a  -y^i  c2> 
the  altitude  interval 


Cg,  Ug  and  iig  ^represent  the  sound  velocity  at  the  lower  and  upper  limits  of 


C-  =  mean  temperature  through  4y1 
Cg  =  mean  temperature  through  AjTg 

Although  Table  m  indicates  that  a  negative  temperature  gradient  exists  in  Ay^  if>?  the  atmospheric  models 
of  example  (b),  the  positive  temperature  gradient  which  has  been  postulated  as  a  necessary  condition  for 
.destructive  focusing  appears  in  the  altitude  interval  4y^. 

3)  Example  (c)  Figure  7c. 

A  dry  adiabatic  temperature  gradient  was  assumed  to  extend  from  the  ground  to  1000  ft.  (-5.4°F/1000 
ft.)  Further,  a  wind  velocity  gradient  of  +  10  ft/sec.,  with  an  initial  velocity  m  =  0,  was  assumed  for  the  same 
altitude  interval.  It  was  found  that  a  sound  ray  of  Initial  angle  0  =  5.232°  would  be  just  refracted  in  the  1000 
ft.  interval.  Therefore  sound  rays  of  initial  angles  0  =  6°  and  0  =  7°  were  investigated  to  determine  the 
atmospheric  conditions  above  1000  ft.  necessary  to  destructively  focus  these  rays  with  ray  0  =  5.232°. 

As  might  be  expected  from  the  results  already  discussed  in  the  two  preceding  examples,  the  mean 
temperature  In  altitude  interval  Yg  for  rays  0  =  6°  and  0  =  7°  had  to  be  greater  than  the  mean  temperature  in 
A  y1  and  this  specified  a  positive  temperature  gradient  of  no  mean  proportions  in  the  altitude  interval  Ay% 
(dependent  upon  the  sound  i  ..y  Investigated  and  greater  than  5,4°  F/1000  ft,  for  0  =  6°).  Thus  it  is  indicated 
that  in  the  normal  type  atmospheric  structure  where  the  temperature  decreases  with  altitude,  a  considerable 
temperature  inversion  is  necessary  to  produce  destructive  focusing  of  sound  rays. 


